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ABSTRACT 

We present color-magnitude and morphological analysis of 54 low-redshift ultralu- 
minous infrared galaxies (ULIRGs; 0.018 < z < 0.265 with z me ^ ian = 0.151), a subset of 
the IRAS Uy sample (Kim & Sanders, 1998), in the Sloan Digital Sky Survey (SDSS). 
The ULIRGs are both bright and blue: they are on average 1 magnitude brighter in 
Mo.i r than the SDSS galaxies within the same redshift range, and 0.2 magnitude bluer 
in 01 (/ — 01 r. They form a group in the color-magnitude diagram distinct from both 
the red sequence and the blue cloud formed by the SDSS galaxies: 24 out of the 52 
unsaturated objects (~ 46%) lie outside the 90% level number density contour of the 
SDSS galaxies. The majority (47, or ~ 87%) have the colors typical of the blue cloud, 
and only 4 (~ 6%) sources are located in the red sequence. While ULIRGs are popularly 
thought to be precursors to a QSO phase, we find few (3 or ~ 6%) in the "green valley" 
where the majority of the X-ray and IR selected AGNs are found. Moreover, none of the 
AGN-host ULIRGs are found in the green valley. For the 14 previously spectroscopic 
identified AGNs (~ 28%), we perform PSF subtractions and find that on average the 
central point sources contribute less than one third to the total luminosity, and that 
their high optical luminosities and overall blue colors are apparently the result of star 
formation activity of the host galaxies. Visual inspection of the SDSS images reveals a 
wide range of morphologies including many close pairs, tidal tails, and otherwise dis- 
turbed profiles, in strong support of previous studies and the general view of ULIRGs as 
major mergers of gas-rich disk galaxies. A detailed morphology analysis using Gini and 
M20 coefficients shows that slightly less than one half (~ 42% in g band) of the ULIRGs 
are located in the merger region defined by morphology studies of local galaxies, while 
the remaining sources are located in the region of late-type and irregular galaxies. The 
heterogeneous distribution of ULIRGs in the G — M20 space is qualitatively consistent 
with the results found by numerical simulations of disk-disk mergers, and our study 
also shows that the measured morphological parameters are systematically affected by 
the signal-to-noise ratio and thus the merging galaxies can appear in various regions 
of the G — M20 parameter space. We briefly discuss the origins of the uncertainties 
and note that the morphology measurements should be implemented with caution for 
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low physical resolution images. In general, our results reinforce the view that ULIRGs 
contain young stellar populations and are mergers in progress, but we do not observe 
the concentration of ULIRGs/ AGN in the green valley as found by other studies. Our 
study provides a uniform comparison sample for studying dusty starbursts at higher 
redshifts such as Spitzer MIPS 24//m-selected ULIRGs at z = 1 ~ 2 or submillimeter 
galaxies. 

Subject headings: galaxies:evolution - galaxiesdundamental parameters - galaxies interactions 
- galaxies:Seyfert - galaxies:starburst - infrared: galaxies 

1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs; Ljr > 1O 12 L0) were first discovered in a large num- 
ber by Infrared Astronomical Satellite (IRAS) two decades ago (Rowan-Robinson et al. 1984; Soifer 
et al. 1984). These objects are among the most luminous sources in the universe, with most of 
their energy radiated in the infrared. Although the powering sources of ULIRGs remain uncertain, 
there is evidence that both AGNs and dusty starbursts contribute to their high bolometric lumi- 
nosities (see Sanders & Mirabel 1996; Lonsdale, Farrah, & Smith 2006 for reviews). The majority 
of local ULIRGs are mergers with tidal features commonly seen (e.g. Sanders et al. 1988; Melnick 
& Mirabel 1990; Murphy et al. 1996), and distant ULIRGs are the possible progenitors of today's 
massive ellipticals (e.g. Barnes & Hernquist 1996). 

Recent large optical photometric surveys show that the color distribution of galaxies are bi- 
modal (Strateva et al. 2001; Hogg et al. 2003; Blanton et al. 2003c; Bell et al. 2004a,b). 
Color-magnitude diagram reveals a red sequence of early type galaxies with little star formation 
activity, and a blue cloud of late type galaxies that form stars. The characteristic number density in 
the standard luminosity function of Schechter (1976), 0*, has doubled for the red sequence galaxies 
since z ~ 1, but hardly changes for the blue cloud, implying that the number and stellar mass of 
the red sequence galaxies have been gradually built up, whereas those of the blue cloud galaxies 
remain nearly constant (e.g., Bell et al. 2004b; Faber et al. 2007). One scenario of building up 
the most massive galaxies in the red sequence is major merging of disk galaxies (Bell et al. 2004b), 
and thus galaxies in the "green valley", the region between the red sequence and blue cloud in 
the color-magnitude diagram, may represent this transitional type. As ULIRGs are proposed to 
be disk-disk mergers that will finally evolve to massive ellipticals, studying the color-magnitude 
relation of ULIRGs may shed light on the evolutionary scenarios of building up massive ellipticals. 

Galaxy morphologies indicate the distribution of baryonic matter in the galaxies, and are 
therefore tracers of its formation and evolution. The majority of local ULIRGs are merging sys- 
tems in which disturbed morphologies such as tidal features are commonly observed. Recent deep 
sky surveys have generated enormous amounts of galaxy imaging data, and automated morphology 
measurement methods are required to analyze those large datasets with much higher efficiencies 
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than visual inspections provide. There are two quantitative approaches to describing galaxy mor- 
phology: parametric methods and non-parametric methods. Parametric methods fit the light dis- 
tribution of a galaxy to pre-defined formulae, such as Sersic index fitting (Blanton et al. 2003c) and 
bulge-to-disk ratio (Peng et al. 2002; Simard et al. 2002), and are usually inadequate to describe 
the morphologies of complex systems. Non-parametric methods, such as the concentration, asym- 
metry, and dumpiness (CAS) system (e.g. Isserstedt k, Schindler 1986; Abraham et al. 1994; Wu 
et al. 1999; Conselice et al. 2000; Conselice 2003), do not rely on pre-defined functions, and may 
be the better methods for describing complex systems. Abraham et al. (2003) introduced the Gini 
coefficient (G), which describes the relative intensity distributions of a galaxy. The Gini coefficient 
is correlated with concentration and increases with the fraction of light in compact components. 
Lotz, Primack, & Madau (2004; hereafter LPM04) introduced M20, which is the second moment 
of the brightest 20% of galaxy flux. M20 is sensitive to spatial distributions of off-axis clumps. G 
and M20 thus should be sensitive to the presence of merger features and indeed LPM04 found that 
local ULIRGs are well separated from normal galaxies in G — M20 space, in the sense that ULIRGs 
have higher G and M20 values than normal galaxies. 

In this paper we present the results of color-magnitude analysis and G and M20 morphological 
analysis of 54 ULIRGs from the IRAS Uy ULIRG sample (Kim 1995; Kim & Sanders, 1998) 
that were observed by Data Release 5 of Sloan Digital Sky Survey. The goals of this paper arc 
to place ULIRGs in the context of normal galaxies by comparing their color-magnitude relations 
and morphologies with those of optically-selected SDSS galaxies, and to construct a low-redshift 
comparison sample for studying high-redshift dusty starbursts, such as Spitzer-identified ULIRGs 
and submillimeter galaxies (SMGs). The IRAS Uy ULIRG sample is a complete flux- limited 
sample from the IRAS Faint Source Catalog (Moshir et el. 1990) with S^oum > Uy and 5 > 
—40°, \b\ > 30° (Kim & Sanders 1998). The sample contains 118 sources with a redshift range 
between 0.018 and 0.268 and a median redshift of 0.145. The infrared luminosities range between 
10 12 - 00 and 1O 12 - 9O L , with a median of 10 12 - 19 L Q . They are the most IR luminous galaxies in 
the low redshift universe, and their uniform selection and completeness make them ideal for a 
statistical analysis. Almost all ULIRGs in the sample show visual signs of interactions or mergers, 
but most are at advanced stages, and the median R-band luminosity is ~ 2L* (Kim, Veilleux, & 
Sanders 2002, hereafter KVS02; Veilleux, Kim, & Sanders 2002, hereafter VKS02). Spectroscopic 
study showed that the Seyfert fraction increases with infrared luminosity, and ~ 50% of the galaxies 
with Ijr > 1O 12,3 L present Seyfert characteristics. These AGN-powered ULIRGs have bolometric 
luminosities and near-IR spectra reminiscent of the local quasars. Another 30% of the whole sample 
are classified as Hii-region and ~ 40% are classified as LINERs, in both of which the energy sources 
are thought to be stellar origin from recent starbursts (< few times 10 7 yr) or shocks (Veilleux, 
Kim & Sanders 1999, hereafter VKS99; Kim, Veilleux, and Sanders 1998) Near-IR study using 
2MASS showed that Seyfert galaxies among the Uy sample have much redder colors and steeper 
spectral indices in the near-infrared than the rest, which indicates powerful AGNs as the energy 
sources (Chen & Zhang 2006). Recent 12 CO (J = 1 -»• 0) observations of 17 ULIRGs in the Uy 
Sample, along with 12 other local ULIRGs showed that their large IR luminosity and gas mass 
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to FIR luminosity ratios are consistent with a model where most of their luminosity is powered 
by a brief surge in star formation rate associated with the rapid gas inflow during the merger 
phase, although there is also evidence for powerful AGNs with extreme Lfir/ L' co among a subset 
and they are possibly transitioning to a quasar phase (Chung et al. 2009). In general, previous 
studies suggest an evolutionary sequence of ULIRGs in which ULIRGs start with cool infrared 
colors {hhpm/ fwixm < 0.2), go through a warm ULIRG state (i^m/ZeOjum > 0.2), and finally 
evolve to quasars (Sanders 1988; VKS02). 

In £j2] we present the sample selection. The image analysis is presented in £j3j We present 
the color magnitude relation of the ULIRGs in £JH their morphological analysis in £j5j and the 
summary in Sj6) We use a ACDM cosmology with Hq = 70kms _1 Mpc~ 1 , f2 m = 0.3 and £l\ = 0.7 
throughout the paper, except for the absolute magnitudes, for which we use h = 1, where h = 
-ffo/lOOkms^Mpc -1 , for a direct comparison with existing studies. 

2. SAMPLE SELECTION 

The sample consists of all the sources from the IRAS Uy ULIRG sample (Kim 1995; Kim & 
Sanders 1998) that have been imaged in the fifth release of the Sloan Digital Sky Survey (SDSS 
DR5; Adelman-McCarthy et al. 2007). 

The SDSS is the largest homogeneous, wide-area multiwavelength sky survey available for 
studying galaxy morphology and color-magnitude distributions in the low redshift universe. We 
downloaded SDSS DR5 images at the coordinate position of each source in the Uy ULIRG sample, 
and did a one-by-one comparison with published optical and near-IR images (KVS02). This yielded 
a final sample of 54 ULIRGs identified in SDSS. The 54 sources span a redshift range between 0.018 
to 0.265, with a median of 0.151, and an infrared luminosity range between 10 12 - 00 and 10 12 L 
with a median of 10 12 ' 23 Lq. The mean values of our sample are very close to those in the whole 
Uy sample. We perform Kolmogorov-Smirnov (K-S) tests between our sample and the whole Uy 
sample for the distributions of redshift, infrared luminosity, and infrared colors (i^S/wn / Fm^m) , and 
the probabilities that the two samples are drawn from the same distribution are 80%, 97%, and 
99%, for redshift, infrared luminosity, and infrared colors, respectively. Thus we conclude that our 
sample is a good representative subset of the whole Uy sample. 

The redshifts and the FIR luminosities of the sample are presented in the second and third 
column in Table [1] respectively. There are 14 ULIRGs harboring one or more AGNs according to 
optical and near-IR spectroscopic study (VKS99 and references therein), and they are identified 
accordingly in Table [TJ The AGN sample includes all galaxies with Seyfert activity but excludes all 
LINERs, which are thought to be low-energy AGNs (e.g. Krolik 1998). Earlier studies of the Uy 
ULIRG sample concluded that the AGN fraction increases with IR luminosity (VKS99; VKS02). 

The comparison sample we choose for studying the color-magnitude relation is the 436,762 
galaxies in the DR4 of the NYU Value Added Galaxy Catalog (NYU-VAGC, DR4; Blanton et al. 
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2005) with spectroscopic redshifts in the same redshift range of our ULIRGs. For studying the 
morphological parameters, the comparison sample we selected consists of the galaxies that appear 
on the same SDSS images of the ULIRGs. The size of each image is 13.51' x 8.98', and there are 
1215, 2059, 2019, and 484 comparison galaxies, in the g, r, i, and z band, respectively. We will 
discuss the selection method in more detail in ^5.11 

3. IMAGE ANALYSIS 

We obtained pipeline-reduced clean images (JpC*.fit) from the SDSS DR5 database. We 
adopted the skylevels in the observation auxiliary files {tsField* .fi£) and subtracted them from 
the images. Each image was checked to ensure the soft skyvalue in the header is correct; for in- 
correct sky values, we estimated the sky levels based on the statistics of multiple small empty-sky 
regions on the same image. 

A mosaic of RGB color images of all 54 sources is shown in Fig. [U (Refer to the electronic 
version for the color image). Each image is 100 pixels across, or ~ 40" at SDSS's pixel scale, and 
centered on the source. The color images were prepared using the method described by Lupton 
et al. (2004); g—, r— , and i— band background subtracted images as blue, green, and red images, 
respectively. The images were smoothed with a median window of 2 pixels wide, aligned using 
r— band images, and rotated so that north is up. The relative scales of red, green, and blue 
images were set to 1:1:2, in order to emphasize the short wavelength radiation from young stellar 
populations. The first visual impression is that the majority of the galaxies show disturbed/irregular 
morphologies, as well as strong color gradients and discontinuities. Note that there is a smeared 
three color source on the image of FSC09539+0857. We checked the SDSS Moving Object Catalog 
(MOC, the third release; Ivezic et al. 2002) and found that it is an asteroid with Unique SDSS 
moving object ID sld059. 

The images are arranged according to the g— band G and M20 coefficients of each source (see 
£j5|), so that for each column G increases towards the top, and for each row M20 increases towards 
the left. The coefficients were calculated using the method described by LPM04 (see §5.11) and are 
printed on each image. Although the images are arranged only relatively in the G — M20 parameter 
space, a clear trend is seen: more disturbed and multiple nuclei sources have larger G and M20 and 
are located toward the upper left corner of the mosaic, and the less disturbed sources, with smoother 
and more symmetric morphology, have smaller G and M20 and tend to be located toward the lower 
right corner. Visual inspection of these images also shows a wide range of morphologies including 
many close pairs, tidal tails, and otherwise disturbed profiles, in strong support of previous studies 
and the general view of ULIRGs as major mergers of gas-rich disk galaxies. 

We computed the magnitude within an elliptical aperture whose semi-major axis is 2 Petrosian 
radii (Petrosian 1976) for each source. The Petrosian radius is defined as the radius of a circular 
aperture on which the mean surface brightness is equal to 20% of the average surface brightness 
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contained inside the aperture. Using an aperture as large as 2 Petrosian radii ensures that most low- 
surface-brightness features are included in the photometry, and the SDSS Petrosian magnitudes we 
use to comparison are also obtained within a 2 Petrosian radius aperture for each source (Blanton 
et al. 2001; Yasuda et al. 2001). For each source, we align the five images using the r-band 
image as the reference. Bilinear interpolation is used during the transformation of each image, and 
flux is estimated to be conserved within the 0.05% level. We run SExtractor (Bertin Sz Arnouts 
1996) on the r-band image and obtain the aperture whose semi-major axis is 2 Petrosian radii, 
and perform the photometry for all the five bands using the same aperture. Some sources have 
disconnected multiple components detected by SExtractor, with each component confirmed to be 
a part of the ULIRG by its redshift based on previous studies (KVS02). For these sources, if the 
components do not have overlapping apertures, we adopt the total flux from multiple components 
for photometry; if they have overlapping apertures, we fit the overlapping region with an ellipse and 
measure the flux within the ellipse, then subtract it from the total flux to obtain the corrected flux 
of the entire source. The airmass and the photometric zero points are obtained from the auxiliary 
files {tsField* .fit) and applied when calculating the magnitudes. The SDSS magnitudes were then 
converted to AB magnitudes using the AB-SDSS magnitudes corrections. 

We find that for the ULIRGs the independently measured apertures are not identical in all 
five bands. This is because many ULIRGs are extended and amorphous, and the extensions of low 
surface brightness features can be different at different wavelengths. The mean relative difference of 
the Petrosian radii measured in g and r band, ((rp :9 — rp jr )/rp :r ), is ~ 8%. We perform the g— band 
photometry using the apertures measured in the g band, and find that the measured magnitudes are 
on average 0.03 mag more luminous than the magnitudes measured using the r— band apertures, 
and the standard deviation of the difference is 0.10 mag. Because we use r— band apertures to 
perform the photometry, we adopt this standard deviation as a typical error of the 01 g — °' 1 r color 
for the ULIRGs. 

Our main purpose is to locate the positions of the ULIRGs in the color-magnitude diagram, 
and compare them with the existing survey data. We therefore adopt the k-corrections to SDSS 
bandpasses shifted to z = 0.1, the median SDSS galaxy redshift, employing the kcorrect 
package developed by Blanton et al. (2003a). Each shifted band is denoted by a superscript on the 
left, e.g., 01 g denotes the SDSS g band shifted to z = 0.1, and 01 g — 01 r denotes the k-corrected 
color between two shifted bands, g and 0ml r. Galactic extinction is corrected using the SFD98 
dust map (Schlegel, Finkbeiner & Davis 1998). We derive the absolute magnitude by applying the 
distance modulus to the k-corrected apparent magnitude. We present the photometry results in 
Table 1. 
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4. COLOR-MAGNITUDE RELATION 

4.1. Color and Magnitude of SDSS Comparison Sample 

In Fig. 2 we plot the 0,1 g — 01 r color against k-corrected ai r-band absolute magnitude Mo.i r 
for all comparison sample galaxies, which consists of the 436,762 galaxies from NYU SDSS Value 
Added Catalog (DR4) within the same redshift range. Individual galaxies are shown as small dots 
while the contours represent the galaxy number density. The bimodal distribution of the SDSS 
comparison galaxies is clearly shown in the color-magnitude diagram where the red sequence and 
the blue cloud form two separate groups. We adopt the empirical definition of the red sequence from 
Weinmann et al. (2006), such that the red-sequence galaxies follow the color-magnitude relations, 

0.1 g _o.i r > o.7 - 0.033(Mo.i r - hlogh + 16.5) (1) 

This relation is shown in Fig. 2 where the red sequence lies above the upper straight line. The 
centroids of the red sequence and the blue cloud are separated by ~ 0.3 in 01 g — 0-1 r, and between 
the two populations there is a relatively low density region called the "green valley". We define 
the green valley as a strip below the red sequence with a width of 0.1 in ai g — 01 r, shown as the 
region between the two straight lines in Fig. 2. The fiducial width of 0.1 is chosen to be one third 
of the color difference between the centroids of the red sequence and the blue cloud, and is also 
roughly equal to the color difference between the centroid and the blue limit of the red sequence 
(the upper straight line). The low density region between the red sequence and the blue cloud is 
wider at fainter magnitudes and narrower at brighter magnitudes, and our simple cut may miss 
some green-valley galaxies at fainter magnitudes but include more blue-cloud galaxies at brighter 
magnitudes. However, this will not affect our main results as we discuss later. 



4.2. ULIRGs Are Optically Bright and Blue 

Table [2] lists the statistics of the colors and magnitudes of the ULIRGs and SDSS comparison 
sample and their subsamples. We summarize below those optical characteristics of the ULIRGs 
and their subsamples (AGN and non-AGN) and compare them with those of the SDSS comparison 
sample. 

The first obvious trend is that the ULIRGs are very luminous in the optical. The median 
k-corrected absolute magnitude is Mo.i r = —21.4, or equivalently 2.5 L* given Mo.i r = —20.44 for 
an L* galaxy (Blanton 2003b). This is consistent with the median of R— band absolute magnitudes 
for the whole 1 Jy sample (~ 2L*; KVS02). For comparison the median k-corrected absolute 
magnitude of the SDSS galaxies is Mo.i r = —20.4, which is 1 magnitude fainter than that of the 
ULIRGs. The absolute magnitudes of the ULIRGs range from Mo.i r = -19.78 to Afo.i r = -23.23 
(excluding the two saturated sources), a luminosity range of a factor of 25. Fifty out of the 
54 ULIRGs (~ 93%; including the two saturated sources) are more luminous than the median k- 
corrected r— band absolute magnitude of the SDSS comparison sample. The difference in magnitude 
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distribution between the ULIRGs and the SDSS comparison sample is further illustrated in Fig. 
[3l The ULIRG sample has a much narrower distribution than that of the comparison sample. 
The ULIRGs are on average 1 magnitude brighter than that of the entire comparison sample (see 
panel(a)) as well as the individual subsample of the red sequence, the green valley, and the blue 
cloud (panel (c)). At the highest luminosity bins the ULIRG sample outnumbers all three SDSS 
sub-samples in fraction (Mo.i r — 5logh < —21; panel (c)). 

Another striking trend is that the ULIRGs are very blue, with only a few exceptions. Their 
9 ~ 0,1 r colors span a wide range, from 0.04 to 1.5, and the majority (87%) of the ULIRGs 
have optical colors as blue as those of the blue-cloud galaxies. The median °- 1 g — 01 r color of the 
ULIRGs is 0.58 and very close to that of the blue cloud (0.55). The median ° 4 g — 01 r color of the 
red sequence and green valley are 0.95 and 0.78, respectively. The color distributions of the ULIRGs 
and the SDSS galaxies are further illustrated in Fig. [4^a) and (c). Considering galaxy color as a 
function of galalxy luminosity, we limit the r-band absolute magnitudes of the SDSS comparison 
sample to the same magnitude range as that of the ULIRGs (—19.8 < Mo.i r — 5logh < —23.2). 
The peak of the ULIRG distribution is ~ 0.3 mag bluer than that of the SDSS sample within the 
same r-band magnitude range, but very similar to that of the blue cloud within the same r-band 
magnitude range (° g — r ~ 0.6). The ULIRGs have a much larger fraction in the bluer bins 
(0.5 < 0-1 g — r < 0.7) than the comparison sample, but a much lower fraction in the redder 
bins (0.7 < 01 g — °' 1 r < 1.0), except for the reddest bins — °' 1 r > 1.2) and the bluest bins 
( 01 g — 01 r < 0.1). The ULIRG sample outnumbers the blue-cloud galaxies in fraction at blue 
colors (°' 1 (? — r < 0.7), has a similar fraction compared to the green- valley galaxies at green 
colors (0.7 < 01 g — 01 r < 0.9), and has a much lower fraction than the red-sequence galaxies at 
red colors (° g — °' 1 r > 0.9). 

The blue colors of the ULIRGs seem inconsistent with the working hypothesis that ULIRGs 
are dusty, massive systems. However, the patchy mix of blue and red regions in the optical color 
images (see color version of Fig. 1) suggests that the dust extinction is patchy, and the extended 
distribution of young, blue starts dominate the overall color of these galaxies. Veilleux, Sanders &: 
Kim (1999) have discussed the blue optical continuum colors of the ULIRGs and suggested that 
the stellar light suffer less dust extinction than the emission-line gas. Our ULIRGs would appear 
even more luminous and bluer in the optical if corrected for the internal extinction. 

The AGN ULIRGs are among the most luminous sources in the ULIRG sample, with a median 
Mo.i r = —21.9, 0.8 magnitude more luminous than the median of the non-AGN ULIRGs. As shown 
in Fig. E^b) the AGN ULIRGs are more concentrated toward high luminosities and outnumber the 
non-AGN ULIRGs at the highest luminosity bin (Mo.i r — 5logh < —23). They also have a much 
flatter color distribution than the non-AGN ULIRGs (see Fig. IHb)). There is little overlap between 
the AGN ULIRGs and the green- valley galaxies in both magnitude distribution (see Fig. E](d)) and 
color distribution (see Fig. IU(d)). 
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4.3. ULIRGs Are Scarce in the Green Valley 

We overlay the colors and magnitudes of the ULIRG sample on top of the SDSS contours in 
Fig. 2, with different symbols representing ULIRGs known to host an AGN (AGN ULIRGs, circles) 
and Hii-region like and LINER galaxies (non-AGN ULIRGs, stars). The images of FSC12540+5708 
and FSC12265+0219 are saturated in the g and r band, and their magnitudes are shown as upper 
limits. The error bar shows a(gp, g —gp )r ), the standard deviation of difference in g— band magnitude 
measured between g— and r— band apertures ($3]). 

The ULIRGs form a distinct group in the color-magnitude diagram: 24 out of the 52 unsatu- 
rated ULIRGs (~ 46%) lie outside the 90% level contour of the SDSS galaxies in the same redshift 
range. Strikingly, only 3 (6^%3) lie in the green valley, and none hosts an AGN. The majority of 
the ULIRGs (47 out of 54, or 87± 7 n %) are located below the bottom straight line that defines the 
red limit of the blue cloud, and 4 (Tt 6 %) are located above the upper straight line that defines the 
blue limit of the red sequence. 

In the local universe most ULIRGs are major mergers of late-type galaxies, and are proposed 
to be the precursors of massive ellipticals (e.g. Toomre & Toomre 1972; Toomre 1977; Mihos & 
Hernquist 1994, 1996; Barnes & Hernquist 1996). Given the favorable scenario in which the red 
sequence is built up by "wet mergers" of late-type galaxies (e.g. Bell et al. 2004, Faber et al. 
2007), one may expect a concentration of merging galaxies in the green valley, the area connecting 
the blue cloud and the red sequence in the color-magnitude diagram. The low fraction (6%) of our 
ULIRGs in the green valley, however, suggests that most of these ULIRGs are not in the transition 
phase from the blue cloud to the red sequence. Instead, the majority (87%) of these ULIRGs have 
typical blue-cloud colors, suggesting that they are still undergoing rapid star formation, and the 
significant quenching of star formation has yet to start. 

AGNs may play a crucial role in quenching star formation and accelerate the transition of 
blue-cloud galaxies to the red sequence (e.g. Hopkins et al. 2006; Croton et al. 2006). Indeed, 
some studies have found an excess of X-ray selected AGNs in the green valley (e.g. Georgakakis 
et al. 2008; Silverman et al. 2008; Hickox et al. 2009; Schawinski et al. 2009). None of the AGN 
ULIRGs in our sample, however, is located in the green valley, and only two (17%) are redder than 
the blue limit of the red sequence. The remaining 10 (83%) unsaturated sources are bluer than the 
red limit of the blue cloud. 

Hickox et al. (2009) suggested an evolutionary sequence for massive galaxies with halo masses 
between 10 12 ~ 10 13 Mq, in which a high accretion rate optical- and infrared-bright SMG/ULIRG/quasar 
phase evolves quickly (< 10 8 yr), and is followed by a slightly lower accretion rate "green- valley" 
phase (lasting ~ 1 Gyr) when the objects are detected as X-ray AGNs. This phase is followed 



The uncertainty on the fraction of sources in any given region of two-dimensional parameter spaces, e.g. the 
fraction of red-sequence galaxies in the color-magnitude relations, or the fraction of mergers in the G — M20 plot, is 
calculated from the possible miscounted number of sources in this region given typical errors in both parameters. 
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by an intermittent radio AGN phase in which the galaxies reside in the red sequence with lowest 
accretion rates. In the context of the Hickox et al. study, the majority of our AGN ULIRGs are 
still in the early optical-bright phase and have yet to evolve to the green valley. 

4.4. Color and Magnitude of AGN Host Galaxies 

In the optical bands the AGN ULIRGs are among the most luminous in the ULIRG sample, 
and most have blue colors (see §3J). AGNs are point-like sources residing in the nucleus region and 
typically have blue optical colors. Thus one may expect the presence of an AGN would affect the 
color and magnitude of the host galaxy (e.g. Hickox et al. 2009). However our AGN ULIRGs 
are low-luminosity Seyfert galaxies and the AGN emission at optical wavelengths may be highly 
attenuated. In order to test whether their high luminosity and blue color is due to a central AGN 
we estimate the AGN color and luminosity contribution by a simple model fitting. 

4-4- 1- Methodology 

Because the resolution of the images is limited to only ~ 1.3" (FWHM, a physical size of 
~ 2.4 kpc at z ~ 0.1, much larger than the physical size of the nuclear region, we do not try 
to decompose precisely the central AGN and the host galaxy light distribution from each image. 
Instead we subtract a scaled PSF from the nucleus of each galaxy and measure the flux of the scaled 
PSF and compare it with the flux of the host galaxy. As there are few unsaturated field stars on 
each image, we adopt the PSFs from the SDSS archive. We use two different subtraction methods 
and obtain a range of the flux for each subtracted scaled PSF and the host. In the maximum 
subtraction method, the maximally scaled PSF is subtracted from the galaxy image such that the 
residual is minimized but does not have negative pixels within a PSF FWHM size aperture around 
the center. By using this method we obtain an upper limit to the point source contribution and a 
lower limit of the host flux. In the second approach (smooth host subtraction method), the PSF is 
scaled and subtracted such that the host has a smooth profile at the center of the subtraction. The 
smooth host subtraction is based on the assumption that the galaxy is composed of a central point 
source and a smooth host galaxy (e.g. McLeod & McLeod 2001). We choose a 3 by 3 pixels area 
centered on the subtraction center and determine the host to be smooth if the standard deviation 
of these nine pixels is less than the average background rms noise on the same image. For objects 
with multiple nuclei we subtract the scaled PSF from only the brightest nucleus. The exceptions 
are FSC13451+1232 (in g and r band) and FSC15001+1433 (in g band), where the two nuclei are 
less than 4" apart with flux difference less than 20%. In these cases we subtracted scaled PSFs 
from both nuclei. We note that FSC12265+0219 and FSC12540+5708 are saturated in both g- 
and r— band images and are excluded from our analysis. 

The flux of the subtracted point source is sensitive to its center position at the sub-pixel 
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level. To find the center we first fit a two-dimensional gaussian profile for each source, and if the 
fitted FWHM is less than 30% larger than the FWHM of the PSF, we adopt the fitted center as 
the subtraction center. Otherwise we adopt the centroid position as the subtraction center. The 
centroid is defined as the location where the spatial derivative of the pixel intensity becomes zero. 
The photometric uncertainty of the subtracted point source is calculated assuming it is dominated 
by the random noise within the central FWHM size aperture. The mean photometric uncertainties 
at g and r band are ~ 0.05 mag and ~ 0.02 mag for the maximum subtraction method, and 
~ 0.05 mag and ~ 0.04 mag for the smooth host subtraction method, respectively. 



4-4-%- AGN Host Galaxies on the Color-magnitude Diagram 

On average the subtracted point sources contribute only a small amount to the total luminosity 
of the AGN ULIRGs in both g and r band. This suggests that most of their optical luminosity is 
associated with the stellar hosts. The maximum, median, and mean ratio of the subtracted point 
source to the total flux, using the maximum subtraction method, are 67%, 28%, 31% in g band, 
and 74%, 27%, and 31% in r band, respectively. Using the smooth host subtraction method, these 
values are 61%, 11%, and 21% in g band, 74%, 8%, and 21% in r band, respectively. The maximum, 
median, and mean magnitude difference between the hosts and the un-subtracted sources are 1.19, 
0.36, and 0.45 (1.48, 0.34, and 0.47) in g band (r band), using the maximum subtraction method, 
and the corresponding values are 1.02, 0.12, and 0.31 (1.03, 0.09, and 0.31) in g band (r band), 
respectively, using the smooth host subtraction method. 

We overlay the 0A g — r vs. Mo.i r color-magnitude diagram for the 12 AGN ULIRGs (filled 
circles) and their AGN-subtracted hosts (open circles) on the SDSS color-magnitude contours in 
Fig. Some hosts are much bluer after the subtraction (e.g. FSC11119+3257, bluer by 0.42 using 
both methods), some are much redder (e.g. FSC01572+0009, redder by 0.33 using the maximum 
subtraction method). However, the median color difference between the host and the original source 
is only 0.007 and 0.005 mag, using the maximum subtraction method and the smooth host method, 
respectively. Even after the maximum subtraction, only one (FSC11119+3257) source moved closer 
to the green valley, and the remaining 11 sources are located close to their original positions in 
the color-magnitude diagram and barely overlapped with the SDSS comparison galaxies. This 
distinction is clearer using the smooth host subtraction method, after which the hosts are much 
closer to their un-subtracted counterparts. In either case most hosts (10 out of 12) remain luminous 
and blue in the color-magnitude diagram. With or without the point source subtraction, these AGN 
host ULIRGs do not appear in the green valley, the hypothesized transitional stage between the 
blue cloud and the red sequence that was associated with X-ray and IR-selected AGNs (e.g. Hickox 
et al. 2009). The blue colors of our AGN host galaxies are broadly consistent with other studies 
of optically-selected AGNs (e.g. Kauffmann et al. 2003; Jahnke et al. 2004; Silverman et al. 
2008). Kaviraj (2008) studied a local LIRG (L > 10 11 L Q ) sample and also found that virtually 
all LIRGs are in the blue cloud and the AGN has no impact on the star formation in their host 
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galaxies. In a broad agreement with out findings, Cowie & Barger (2008) also found that the bulk 
of extinction-corrected 24^m selected sources in the GOODS-North field lie in the blue cloud. 

Our study is limited by the small sample size, and the point source decomposition is limited 
by the spatial resolution of SDSS. A more complete sample observed with HST or the next gener- 
ation James Webb Space Telescope (JWST; Gardner et al. 2006 ) is necessary in order to fully 
characterize the optical properties of the host galaxies in AGN ULIRGs. 

4.5. Color-magnitude Relation: ULIRGs vs Low-z QSOs and Low-z Type 2 Quasars 

Motivated by the high infrared and bolometric luminosity, it has been proposed that ULIRGs 
are the early stage of dusty quasars, and a ULIRG-QSO evolutionary scenario has been suggested 
(Sanders et al. 1988). Although the sources of the bolometric luminosity for ULIRGs are nearly 
completely obscured by dust in the visible bands, a comparison of their host properties with those 
of QSOs should provide an important test for the proposed evolutionary scenario. For example, if 
the ULIRGs and QSOs are identical except for a geometrical difference such as the viewing angle 
and dust geometry, then the optical properties of their hosts should be largely identical. If a purely 
temporal evolution and an AGN feedback process removing gas and dust are the key difference 
(e.g. Narayanan et al. 2009), then a significant overlap with a systematic shift is expected in the 
host color-magnitude relation. 

The QSO sample we selected is composed of 1070 objects in the SDSS DR5 quasar catalog 
(Schneider et al. 2007) within the same redshift range (0.018 < z < 0.265) as the ULIRGs. The 
SDSS QSOs are type 1 quasars, selected based on their broadband colors, with M, < —22 mag 
and spectroscopically confirmed to have broad lines (FWHM greater than 1000 kms -1 ). We also 
selected a type 2 quasar sample for comparison, which consists of 402 objects in the SDSS type 2 
quasar catalog (Reyes et al. 2008) within the same redshift range, selected based on their optical 
emission lines. Systematic photometric difference should be minimal among the SDSS QSOs, type 
2 quasars, and the ULIRGs due to the shared origin of the data. For SDSS QSOs and type 2 
quasars, we apply the k-correction and extinction correction using the same procedure as for the 
ULIRGs to obtain the absolute magnitudes. We plot in Fig. [6] the 01 g — r colors vs Mo.i r of 
the QSOs and the type 2 quasars with dots and contours, on top of the contours of the SDSS field 
galaxies and the ULIRG symbols. 

Strikingly, the distribution of the type 2 quasars peaks in the green valley where the X-ray 
and IR-selected AGNs are found (e.g. Hickox et al. 2009). In comparison, the SDSS QSOs only 
overlap slightly with the blue cloud and extend toward very blue colors. The difference in their 
distributions in the color-magnitude diagram is at least partly due to selection effects. The SDSS 
QSOs are selected unobscured type 1 quasars where the blue continuum and high luminosity of 
the central AGNs are visible, while the type 2 quasars are selected obscured AGNs and thus the 
dust and gas block the view to the central blue and luminous AGNs. Studying whether these green 
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type 2 quasars represent galaxies in a transitional stage from the blue cloud to the red sequence 
by AGN quenching, or their green colors are simply color combinations of redden obscured AGNs 
and the host galaxies, is beyond the scope of this paper. 

The ULIRGs are distributed very differently from the SDSS QSOs and the type 2 quasars in 
the color-magnitude diagram. The mean absolute magnitude of the ULIRGs ((Mo.i r .) = —21.4) is 
very close to that of the SDSS QSOs ((Mo.i r ) = -21.3), but their mean color ({ 0A g - ai r) = 0.57) 
is 0.35 mag redder than that of the SDSS QSOs (0.22). The SDSS QSOs occupy a narrow range in 
Mo.i r , overlap only slightly with the ULIRGs in the color magnitude diagram, and extend farther 
to bluer colors than the ULIRGs. In comparison, the (Mo.i r ) of ULIRGs is 1.1 mag more luminous 
than that of the type 2 quasar (-20.3), and their ( 0,1 g — 0,1 r) is 0.14 mag bluer (0.57 compared to 
0.74). The ULIRGs overlap with the type 2 quasars only at faint magnitudes (Mo.i r > —21), and 
all but one AGN ULIRGs are more luminous than 90% of the type 2 quasars. We also implement 
two-dimensional Kolmogorov-Smirnov tests (KS test, Peacock et al. 1983; Fasano & Franceschini 
1987) between different samples in the color-magnitude two dimensional parameter space, and the 
results are listed in Table El The test results indicate that neither the ULIRG sample nor the AGN- 
ULIRG subsample is drawn from the same underlying distribution with the SDSS QSO sample or 
with the type 2 quasar sample. 



4.6. Evolutionary Tracks on the Color- magnitude Diagram for ULIRGs 

We model the color-magnitude evolutionary tracks for the ULIRGs using the stellar synthesis 
model BC03 (Bruzual & Chariot 2003) in order to gain some insight into the distribution and 
evolution scenarios of ULIRGs on the color-magnitude diagram. Rather than providing a full 
solution to the SEDs with best-fit parameters, the primary goal of our models is to demonstrate 
qualitatively the possible evolutionary paths these galaxies may follow on the color-magnitude 
diagram and to establish the associated timescales. 

We start with the simplest cases in which the star formation histories are: 1) a passively 
evolved single stellar population burst (SSP), and 2) a constant star formation rate (constant). 
Fig. [7] shows the evolutionary tracks for these two models with a solid line (SSP) and a dashed 
line (constant). Both tracks are k-corrected to z = 0.1 in the same manner as the ULIRGs. The 
SSP track is normalized so that it passes through the median color ( 01 g — r = 0.58) at the 
median absolute magnitude of the ULIRGs (Mo.i r = —21.43, or L r = 2.5L r ^). The constant track 
is normalized to reach 2.5L* at t = 20Gyr. We use solar metallicity in our models for simplicity and 
ignore intrinsic extinction as we are primarily concerned with the unobscured stellar population. 
Three symbols are plotted on each track to mark three ages, 1 Gyr, 2 Gyr, and 5 Gyr. In the SSP 
model, we find that the stellar populations of ULIRGs have ages spanning from 200 Myr to 8 Gyr 
based on their optical colors, with a median age of ~ 1 Gyr. One exception is the extremely red 
source FSC11119-I-3257 whose color is too red to be reproduced by the SSP model. Our modeling 
suggests that a galaxy with a median age takes 1.5 Gyr to reach the red sequence. In the constant 
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star formation case, the reddest color of the stellar population (at t = 20Gyr) is 01 g — °' 1 r = 0.43, 
which is bluer than 80% of the ULIRGs. This implies that in the constant star formation case the 
star formation has to be quenched in order to move the ULIRGs to the red sequence. 

Next we construct a model with a more complex star formation history and apply different 
metallicity and dust treatments to try to understand their effects on the color and magnitude. 
Our model assumes that: (a) a ULIRG is a merger system consisting of two equal-mass late type 
galaxies and that the merger triggers a starburst; (b) the star formation history consists of an SSP 
component associated with each progenitor and an exponentially decaying starburst triggered by 
the merger; (c) the SSP components are identical with the same age and the starburst is triggered 
when the SSP components evolve to the median color and magnitude of the ULIRGs; (d) the 
starburst forms as much as 10% of its progenitors' stellar mass with an e-folding time of 10 7 yr. 
The timescale of 10 7 yr is close to the dynamical timescale of ULIRGs such as Arp 220 (6 x 10 6 
yr; Mauersberger et al. 1996) and 10% is a fiducial fraction of mass formed during a merger that 
is consistent with the results in the numerical simulations of disk mergers (e.g. Mihos & Hernquist 
1996). Given that the progenitors are two L* galaxies with M* ~ 1O 11 M0, the calculated star 
formation rate starts with a maximum value of ~ 2000MQyr~ 1 and decreases to several hundred 
M yr _1 after an e-folding time. Although these star formation rates seem extremely high, they are 
required by the observed high FIR luminosities. We calculate the star formation rate of the ULIRGs 
from their FIR luminosities using the FIR-SFR relation (Kennicutt 1998). The star formation rates 
range from 170 M Q yr _1 to 1370 MQyr -1 , with a median of 270 MQyr -1 , which are consistent with 
our model predicted values. 

We examine three different metallicities, Z = 0.008, Z = 0.02, and Z = 0.05, corresponding to 
a sub-solar, solar, and super-solar metallicity, respectively, and construct three base models using 
this star formation history. We assume that the starburst has the same metallicity as the passively 
evolved progenitor populations. We set up a grid of different dust extinctions following Chariot 
& Fall (2000) for each metallicity, but apply the dust extinction only to the starburst population. 
There are two parameters in the dust treatment method: the total effective V-band optical depth 
Ty that affects stars younger than 10 7 yr, and the fraction \i of extinction that comes from the 
diffuse ISM and thus affects stars of all ages. We choose Ty between 1 to 10 with a step size of 1. 
Bruzual & Chariot (2003) noted that the average value of [i is 0.3, and we adopt this value for our 
modeling. Our choice of optical depth is conservative, based on the common concept that these 
galaxies are dusty starbursts and/or AGNs. However we reiterate that the ULIRGs on average 
are blue and bright, and in order to produce the blue colors the dust attenuation associated with 
the escaping optical emission has to be modest. As discussed earlier, it is also possible that the 
geometry of the dust distribution in the system is patchy so that the high infrared luminosity 
originates from the most dust-attenuated regions, while the optical emission comes from relatively 
less attenuated regions, e.g. tidal structures at large distances. 

We plot nine modeled tracks along with the ULIRGs on the color-magnitude diagram in Fig. 
The tracks start from the time when the starburst is induced. The three rows of panels show 
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sub-solar, solar, and super-solar metallicity, from top to bottom, respectively. For each metallicity 
there are three tracks with ry = 1, 5, 9, from left to right, respectively. The tracks are k-corrected 
to z = 0.1 and normalized to the median Mo.i r which is roughly 2.5 L*. Symbols are plotted on the 
tracks to mark the age of 1, 2, and 5 Gyr after the starburst starts. These tracks clearly suggest 
that both metallicity and dust extinction play important roles in our understanding of the evolution 
of these systems. Dust extinction mostly affects the blue optical colors of the systems younger than 
several hundred million years. For all metallicities, the stellar populations reach their bluest colors 
in less than 20 Myr, but the bluest color a stellar population can reach becomes increasingly redder 
with increasing dust extinction. The blue colors of many ULIRGs can be produced only with 
the lowest dust extinctions. However, there is almost no color difference between different dust 
extinctions when the populations are older than 1 Gyr and their 

o.i g _o.i r 

colors are redder than 

0.7, and any systems with the same metallicity will redden to the red sequence at roughly the same 
time. Metallicity affects the range of colors and the timescale of the color evolution of the stellar 
populations, and higher metallicity systems span wider ranges in the color space, and evolve more 
quickly. As shown in Fig. EJ the bluest ULIRGs can be reproduced only by the highest metallicity 
track. After 20 Gyr, the super-solar track is 0.25 mag redder than the sub-solar case. It takes a 
sub-solar population 2 Gyr to evolve to the red sequence, 1.5 Gyr for the solar population, and 
only 1 Gyr for the super-solar one. 

In reality the star formation history is almost certainly more complicated than the simple 
models we discuss here. In the third model we use a simplified version of star formation his- 
tory produced from numerical simulations of gaseous disk mergers by Springel, Di Matteo, and 
Hernquist (2005). Springel et al. used GADGET-2 (Springel 2005) to trace star formation in a 
merger of two spiral galaxies and studied the effects of back hole feedback on the quenching of 
star formation. The two spiral galaxies have equal dynamical masses of 3.85 x 1O 12 M , and form 
stars at an exponentially decaying rate before the merger, with the pre-merger star formation rate 
~ 200MQyr _1 . The merger-triggered starburst happens at t = 1.5Gyr and the peak star formation 
rate is ~ 2OOOM yr~ 1 . In the black hole feedback scenario the star formation is quenched right 
after the starburst and the star formation rate decreases to less than 1 M^yr -1 at t = 2.5Gyr, 
while in the no feedback scenario the star formation rate decays slowly and remains at several solar 
masses per year for several Gyr. The intrinsic extinction is ignored in this model. We plot the 
evolutionary tracks in Fig. [7]with the dotted line representing the track without AGN feedback and 
the dashed line representing the one with AGN feedback. Three symbols are plotted on each track 
to mark the ages of 1, 2, and 5 Gyr. We find that some of our most luminous ULIRGs are located 
very close to these tracks in the color-magnitude diagram. Both tracks evolve more slowly than 
the SSP track on the same plot because the starburst produces new young blue stars. Between the 
two tracks, the one with AGN feedback clearly reddens more quickly than the one without AGN 
feedback: it takes 8.5 Gyr (or 7 Gyr after the starburst) for the galaxies without AGN feedback 
to evolve to the red sequence, which is 1.5 Gyr longer than the one with AGN feedback. The 
difference in the timescales is consistent with the scenario in which AGN feedback quenches the 
star formation effectively, and consequently helps form the color bimodality of galaxies. 
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In summary, these simple models suggest possible evolutionary paths of the ULIRGs on the 
color-magnitude diagram. We conclude that constant star forming galaxies cannot evolve to the 
green valley without the star formation quenched. Galaxies with lower star formation activity 
evolve to the red sequence more quickly and spend shorter time in the green valley, as seen in the 
SSP model and the numerical simulation model. Dust affects mostly the blue colors of the stellar 
population, and has almost no influences on the red colors and the timescale of evolution to the 
red sequence. Stellar population with higher metallicity explores a wider range in the color space, 
and evolves more quickly to the red sequence. 

The color evolution of the ULIRGs may be presented by any one of the models with different 
masses and dust extinctions. For example, the evolutionary track of the most luminous ULIRGs 
in our sample may be presented by either the numerical simulation models with slightly different 
mass, or by the SSP track in Fig. [7] with ten times more mass; the evolutionary track of a ULIRG 
with the median color and magnitude may be presented by the numerical simulation models with 
90% of the optical emission being attenuated. It is also possible that the color evolution of a galaxy 
is the result of composite star formation histories, so that the galaxy may leave the blue cloud and 
enter the red sequence for more than once. We will not discuss these more complicated scenarios 
because they are out of the scope of this paper. 



5. MORPHOLOGY 

5.1. Gini and M20 Calculations 

We use the code developed by Lotz et al. (LPM04) to calculate the morphology coefficients 
G and M20 for both the 54 ULIRGs and the comparison sample, galaxies appearing on the same 
SDSS images of the ULIRGs. The Gini coefficient, G, is defined as 

n 

G = nm — tt E( 2 * -n-i)\Xi\ (2) 

X n(n — 1 

where n is the number of pixels and X^s are the sorted pixel intensities in increasing order. The 
total second-order moment M to t of the pixels assigned to the galaxy is defined as 

n 
i 

where fi is the flux in each pixel in decreasing order, Xi and are the coordinates of the pixel, and 
x c , y c is the center of the galaxy, computed by minimizing Mtot- M20 is the normalized second-order 
moment of the brightest 20 percent of the galaxy's flux, such that 

M20 = logi (%^) while V ft < 0.2/tot, (4) 
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where jtot is the total flux of the galaxy. 

As noted by LPM04 and Lotz et al. (2008b), the quantity G is very sensitive to the ratio of 
low surface brightness pixels to high surface brightness pixels, so a well-defined segmentation map 
is essential to measure G and M 2 o accurately. Also, as pointed out by Lisker (2008), measurement 
uncertainties of G are minimized when using the Petrosian radius as the aperture size. We thus 
adopt the same approach as LPM04 and assign the pixels brighter than the surface brightness at 
the Petrosian radius (semi- major axis as measured in an elliptical aperture) to the galaxy. The 
measured coefficients are not robust for a source with the signal-to-noise ratio per pixel within the 
Petrosian radius ((S/N)) less than 2.5, or with a Petrosian radius (rp) less than twice the size of 
FWHM of the PSF (~ 1.3" in SDSS) (LPM04 and Lotz et al. 2008b), where the signal-to-noise 
ratio per pixel is defined as the ratio of the mean flux and the mean rms noise within an aperture. 
We therefore exclude all sources with (S/N) < 2.5, and/or rp < 2.6". The u— band images are 
noisy in general: almost 60% of the ULIRGs have (S/N) less than 2.5, and thus we do not include 
the u — band results further in our analysis. We also exclude possible candidates for stars in the 
comparison sample by rejecting sources with the S Extractor star /galaxy separation parameter 
CLASSJ3TAR greater than 0.9. The final selection yields the measured morphology parameters 
for 50, 50, 51, and 42 ULIRGs, and 1215, 2059, 2019, and 484 comparison galaxies, in the g, r, i, 
and z band, respectively. 



5.2. Morphology of the ULIRGs 

In Fig. [9]we plot the g—, r— , i— , and z— band G against M 2 o of the comparison galaxies using 
contours and dots, where the contours are 10 equally spaced levels between 10% and 90% of the 
peak number density, and dots show galaxies located outside the lowest 10% contour level. On 
top of them we overlay the morphology parameters of ULIRGs with circles (AGN ULIRGs) and 
stars (Hii-region like/LINERs). We also plot the empirical lines that divide the parameter space 
into four morphology regions: mergers, elliptical galaxies, irregular galaxies, and spiral galaxies, 
according to the classifications for z ~ galaxies described in LPM04 and Lotz et al. (2008a): 

Merger : G > -0.115M 20 + 0.384 
Elliptical : G < -0.115M 20 + 0.384 k G > O.U5M 20 + 0.769 
Irregular : G > 0.115M 20 + 0.697 k G < 0.115M 20 + 0.769 (5) 
k G < -0.115M 20 + 0.384 
Spiral : G < -0.115M 20 + 0.384 k G < 0.115M 20 + 0.697 

The distributions of the ULIRGs at all wavelengths are consistent in that they form a similarly 
heterogeneous group in this parameter space. The centroids of the distribution in G — M 2 o are 
[0.55,-1.42], [0.55,-1.43], [0.57,-1.47], [0.54,-1.45], in g, r, i, and z band, respectively, all of which 
are located very close (AG < ±0.02 for any given M 2 q) to the separation line between merger 
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and non-merger galaxies. There are only 21(421^%), 21(42±|°%), 25(49^%), and 15(361^%) 
ULIRGs located in the merger region, in g, r, i, and z band, respectively. These are very low 
fractions compared to the z ~ results where the majority (80%) of ULIRGs are located in the 
merger region (LPM04). The distributions of the ULIRGs in the parameter space vary slightly 
with wavelengths: the merger fraction is the highest (49%) using the i— band derived parameters, 
and lowest (36%) using the z— band derived ones. Although the morphology parameters of the 
ULIRGs span a wide range in the parameter space, their G — M20 distribution is clearly different 
from that of the comparison sample galaxies. The distribution centroids of the comparison galaxies 
in the four bands are all located in the empirical region of spiral galaxies, and only 56(5^ 4 %), 
117(6i^ 9 %), 146(7if%), and 26(5+ 21 %) galaxies are found in the empirical region of mergers, in 
g, r, i, and z band, respectively. 

The poor physical resolution of the SDSS images (700 pc pixel -1 at z ~ 0.1) may lead to higher 
random and systematic uncertainties in G and M20 (see LPM04), and thus careful estimation of 
the uncertainties in G and M20 is necessary. Using image and noise simulations (described below in 
§5.3p . we estimate the uncertainties of morphology measurements (AG < 0.05 and M20 < 0.3 — 0.4) 
in SDSS images and plot the largest uncertainties as error bars in Fig. [9j We find that a galaxy seen 
against higher background noise tends to move to the left lower corner in the G — M20 plot, and we 
plot a vector to illustrate the typical changes in G and M20 resulting from adding background noise 
until the source only barely satisfies the signal-to-noise ratio and Petrosian radius criteria. The 
vector shows that a source selected as a merger at high signal-to-noise ratio could move away from 
the merger region and stay in the late-type galaxy region if it is immersed in higher background, and 
consequently, the merger fraction selected by G and M20 is likely to be a lower limit. Although less 
than half of the ULIRGs have empirical merger G — M20 relations, the heterogeneous distributions 
are qualitatively consistent with what was found in the recent numerical simulations of merging 
galaxies of Lotz et al. (2008b). We will discuss these noise simulations further, together with the 
uncertainties in measuring the morphologies, in §5.31 

AGN ULIRGs are slightly more concentrated in the merger region than non-AGN ULIRGs: 
there are 45%, 50%, 60%, 50% of AGN ULIRGs in the merger region, in g, r, i, and z band, 
respectively, compared to 41%, 40%, 48%, 30% for non-AGN ULIRGs. However the difference is 
not statistically significant given the small number of AGN ULIRGs (11 in g band and 10 in other 
bands based on the signal-to-noise ratio and size selection criteria) . There is no significant difference 
between the mean G ((G)) and M 20 ((M 20 » for the AGN ULIRGs and non-AGN ULIRGs: for 
AGN ULIRGs, (G) is 0.58, 0.57, 0.59, 0.58, and (M 20 ) is -1.48, -1.48, -1.36, -1.48, in g, r, i, and 
z band, respectively; in comparison, (G) is 0.54, 0.54, 0.56, 0.53, and (M20) is -1.42, -1.45, -1.49, 
-1.42 for non-AGN ULIRGs in the same bands. 



5.3. Are the ULIRGs Mergers in G - M 20 space? 



Our visual inspection shows that almost all the ULIRGs are disturbed systems (Fig. [2|), and 
this is consistent with the visual classifications by VKS02, who found that all but one of the 118 
IRAS Uy sample show visual merger features. We perform the G — M20 analysis and find that 
more disturbed galaxies do have higher G and M20 coefficients. This is illustrated in Fig. [T] where 
more disturbed galaxies tend to occupy the upper left corner of the mosaic with higher values 
of both G and M20 However our studies also show that ULIRGs are a heterogeneous group in 
G — M20 space: only slightly less than half of the sources lie above the solid line in each panel of 
Fig. [9j the region where most local mergers and ULIRGs were found by LPM04. Interestingly, 
there is only one source located in the early type region, and all the remaining sources fall in the 
region where irregular and late-type galaxies are located. This heterogeneous distribution seems 
to be inconsistent with the result of LPM04 that 80% of local ULIRGs are located in the typical 
merger regions in the parameter space. However, the spread in the parameter space we found is 
supported by numerical simulations. Lotz et al. (2008b) analyzed the morphological parameters in 
SDSS g— band images of mergers of equal mass gas-rich spirals by using the N-body/hydrodynamic 
simulation code GADGET (Springel, Yoshida, & White 2001) and Monte Carlo radiative transfer 
code SUNRISE (Jonsson 2006; Jonsson et al. 2006), and find that the mergers are most disturbed 
in G — M20 at the first pass, but they can have normal galaxies morphologies at other merger stages. 
They also noted that two-thirds of the ULIRGs in LPM04 exhibit double or multiple nuclei, and 
therefore are more effectively selected by G — M<iq- The timescale is only 0.2 — 0.6 Gyr for a 
merger to appear in the empirical merger region defined by LPM04, and the range in this time 
scale depends on dust, orbit parameters, and observing orientations. Thus at face value, our results 
imply that about half of our ULIRGs have been captured within this time window. The fact that 
the other half of our ULIRGs do not appear in the empirical merger region does not mean that 
they are not mergers. Instead, they might be in other merger stages when their morphological 
parameters are consistent with those of normal galaxies. The heterogeneous picture of our ULIRGs 
thus may represent different evolutionary epochs of the ULIRGs. 

In order to gain some insights into the G — M20 morphology we compare our g-band G — M20 
classification with the interaction classification for the same sources made by Veilleux et al. (2002). 
Veilleux et al. classified each object into one of the six sequential merging stages according to 
its morphology. The comparison of two different classifications is illustrated in Fig. [TUJ We find 
that all triple mergers by the classification of Veilleux et al. are also classified as mergers by 
G — M20 Sources in their type Ilia group, i.e., wide binary pre-merger systems with apparent 
separations greater than 10 kpc, are also very likely to be recognized as G — M20 mergers: ~ 80% 
of such sources are recognized as mergers using G — M20 Sources in their other groups are less 
likely to be recognized by G — M20 as mergers, and the fraction of G — M20 mergers appears to 
decrease toward later merging stages. We find that G — M20 mergers make up ~ 31% of the 
close binary pre mergers (type Illb) , ~ 50% of the compact mergers (type IVa), ~ 46% of the 
diffuse mergers (type IVb), and only ~ 11% of the old mergers (type V) are G — M20 mergers, 
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respectively. Sources classified as close binary pre-merger systems (type Illb) and old mergers (type 
V) have very large fractions to be recognized as spiral galaxies by G — M20 (~ 54% and ~ 67%, 
respectively). These comparison results are qualitatively consistent with the simulation results by 
Lotz et al. (2008) in which the most disturbed G — M20 morphology happens during the first pass 
and maximum separation. However, since our measured morphology is limited by relatively small 
number statistics and relatively large uncertainties, we do not conclude strongly that the observed 
morphology can be well explained by the simulations. 

We also find that uncertainties in measuring the morphologies contribute significantly to the 
distribution of the ULIRGs in the parameter space. LPM04 reported that the measurements of 
G and M20 are robust to within 10% when the source has signal-to-noise ratio per pixel {S/N) 
greater than 2.5, the Petrosian radius larger than 2.5 times the PSF FWHM, and the physical 
resolution (parsec per pixel) higher than 500 pc pixel -1 . However when the physical resolution is 
lower the uncertainties of morphologies largely increase (Fig. 6 in their paper) since small structures 
are washed out. At z ~ 0.1 — 0.2 the spatial resolution of the SDSS ULIRGs is roughly 0.7-1.3 
kpc pixel -1 and the uncertainties need to be calibrated carefully. 

We perform simple simulations to estimate the uncertainties of G and M20 for the SDSS 
ULIRGs at different noise levels. We use four model elliptical galaxies, four spiral galaxies selected 
from SDSS images, and three merging galaxies selected from our ULIRG sample. The four elliptical 
galaxies are modeled using de Vacoulers profiles and are added to real SDSS images to mimic the 
observational conditions. The three spiral galaxies are selected with bright g-band magnitudes 
of 13.9, 15.2, and 16.4, respectively, much brighter than the median g-band magnitudes of the 
SDSS sample (18.0 mag). The three ULIRGs (FSC08572+3915, FSC14060+2919, FSC14121-0126) 
are selected with bright g-band magnitudes and unmistakable merger morphology. The g-band 
magnitudes are 16.4, 16.7 and 17.8, respectively, all brighter than the median g-band magnitudes 
of the ULIRGs (17.8 mag). For each image we generate random noise at a series of different levels, 
and at each level we generate 20 noise maps independently and add them to the original image. 
Morphologies, Petrosian radius, and (S/N) are measured for each noise-added source and compared 
to the original measurements. Although our simulation is limited by the sample size, ~ 90% of the 
noise-added sources selected by the signal-to-noise and size criteria ({S/N) > 2.5 and R p > 2.6") 
have g— band magnitudes between 16 and 19, and g— band Petrosian radii between 3 and 12 arcsec 
(~ 6 - 22 kpc at z ~ 0.1), ranges similar to those of the ULIRGs. 

We find that the measured Petrosian radius and the (S/N) do not always decrease monotoni- 
cally with the amount of noise added, and they are broadly distributed at low signal-to-noise levels. 
We also find that, rather than remaining fixed for a given source, the measured G decreases and 
M20 increases systematically with increasing noise. Therefore when seen against higher background 
noise, a galaxy tends to move to the lower left corner on the G — M20 plot. This effect will de- 
crease the fraction of the mergers observed in the merger region. The Gini coefficient G decreases 
by as much as 0.2 with increasing background noise until the galaxy is no longer detected (either 
r p < 2 x FWHM, or (S/N) < 2.5), and M20 increases by as much as 0.5. The standard deviations 
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of G at every noise level are all less than 0.05, and the spiral galaxies have the smallest dispersions. 
The standard deviations of M20 are all less than 0.3-0.4, and the elliptical galaxies have the smallest 
dispersions. There is no systematic difference for the uncertainties between different bands. These 
uncertainties (AG = 0.05 and AM = 0.5) are shown in Fig. [9] as error bars, and the systematic 
trend with increasing noise in G and M20 is shown as an arrow. 

The aperture within which the morphologies are measured is also an important source of un- 
certainty. As noted by Lisker (2008), G value measured within larger apertures have systematically 
larger values, and the measured uncertainties are minimized when the Petrosian radius is used as 
the aperture size. The reason behind is that more low surface brightness pixels are included within 
a larger aperture size, which steepens the intensity distribution of the pixels and consequently in- 
creases the value of G. Some of the low surface brightness features of the source, often at larger 
distances from the center, are mistakenly excluded when the aperture size is chosen too small, which 
flattens the pixel intensity distribution and consequently decreases the value of G. When assigning 
pixels to the sources we adopt the surface brightness at the Petrosian radius as the threshold, and 
assign pixels brighter than the threshold to the source. Therefore, according to their study our 
measured morphologies should have minimized uncertainties, because our apertures should be very 
close to the Petrosian radius and not significantly larger or smaller. 



5.4. Relations among Morphology, Optical Color, and FIR Luminosity 

The morphology, color-magnitude relation, and the infrared luminosity of ULIRGs present a 
rather complicated picture of these systems. The huge infrared luminosity in our ULIRGs suggests 
that their dominant powering sources should be obscured by dust. However, the majority of them 
are optically bright and blue. This further implies that dust is not uniformly distributed and that 
significant amount of the unobscured stellar light is seen directly. Indeed, as shown in Fig. 1, there 
are obvious color gradients in the ULIRGs, with their tidal features at large distances appearing 
blue and their central regions appearing red. Lotz et al. (2008b) suggested that merging systems 
are most disturbed in G — M20 space and are located in the merger region in the G — M20 plot during 
the first pass when the tidal features are prominent. After the first pass they will move to the late- 
type and irregular region (see Fig. 5 in their paper). Thus we would expect the integrated colors 
of those disturbed systems selected by G — M20 to be relatively bluer than the ones with normal 
G — M20 relations, due to the blue colors of the tidal features. Recent hydrodynamic simulations of 
disc-disc mergers also suggested that the infrared luminosity of ULIRGs will peak during the final 
merger when enough metals and dust have been accumulated (e.g. Jonsson et al. 2006), so one 
would also expect sources with higher infrared luminosity are preferentially found in the late-type 
and irregular region. 

We test these hypotheses by plotting the g— band Gini vs M20 for the ULIRGs in Fig. [TT| with 
different symbol sizes representing their FIR luminosities and different symbol colors representing 
their 

o.i 5 _o.i r 

colors. There are no strong relations between the optical colors and G or M20, except 
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that optically redder sources are slightly more preferentially located in the non-merger regions with 
lower G values. This is further illustrated in the small panel in the same plot, where we plot the 
color distributions of sources with G > 0.55, the median G of the whole sample, and of sources 
with G < 0.55. The sources with lower G values tend to distribute toward redder optical colors, 
although the reddest source has a higher-than-median G value (0.62). All AGN ULIRGs have G 
values greater than 0.5. There are no obvious correlations between the FIR luminosity and G or 
M 20 . 

These results are broadly consistent with the scenario that the blue optical colors originate from 
more disturbed systems, and these systems are not during the phase when the FIR emission peaks. 
However, as mentioned earlier, the uncertainties in measuring G and M20 could be substantial, and 
the measured morphologies are affected by orientation, orbital parameters, viewing angle and dust 
content (Lotz et al. 2008b). We thus do not conclude any strong relations among the morphologies, 
the FIR luminosity, and the optical color of these system. 

6. SUMMARY 

We present color-magnitude and morphological analysis of 54 ultraluminous infrared galaxies 
(ULIRGs) in the Sloan Digital Sky Survey (SDSS). The sample consists of all the ULIRGs from the 
IRAS Uy sample (Kim & Sanders, 1998) that have been imaged in the Data Release 5 of SDSS, 
spanning a redshift range from 0.018 to 0.265 with a median redshift of 0.151. The main results 
are: 

• The ULIRGs are a very luminous group of galaxies in the optical, and the majority (~ 93%) 
are more luminous than the median r-band absolute magnitude of the SDSS comparison galaxies 
in the same redshift range. 

• The ULIRG sample forms a distinct group in the color-magnitude diagram. 24 out of the 52 
unsaturated ULIRGs (~ 46%) lie outside the 90% level contour of the SDSS galaxies in the same 
redshift range (NYU-VAGC, DR4; Blanton et al. 2005). The majority of the ULIRGs (~ 87%) 
have typical colors of the blue cloud, only ~ 6% are located in the green valley, and ~ 7% are 
located in the red sequence, which implies that most of the ULIRGs are still undergoing copious 
star formation. 

• There are 14 ULIRGs known to harbor an AGN (AGN ULIRGs). None of them are located 
in the green valley. We use two simple approaches to estimate the point source flux in these 
systems, and find that on average the central point source contributes less than one-third to the 
total luminosity in r band. The host galaxies do not have a significantly different distribution in 
the color-magnitude diagram after the point sources removed. 

• We perform G and M20 analysis on the ULIRGs and find that their distribution in the G—M20 
space is heterogeneous. Less than ~ 50% (e.g. ~ 42% in g band) of the ULIRGs are located in 
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the merger region found for local ULIRGs (LPM04). However the heterogeneous distribution is 
consistent with the morphology produced by numerical simulations (Lotz et al. 2008b), and we 
also discuss the uncertainty in the morphology measurements at the physical resolution of SDSS. 

• We find that the ULIRGs have comparable optical luminosity as the SDSS QSOs within 
the same redshift range but much redder. In comparison, the ULIRGs are much more luminous 
than the SDSS type 2 quasars within the same redshift range and much bluer. The distribution 
of the SDSS type 2 quasars peaks at the green valley. We perform two-dimensional K-S tests and 
the results show that the ULIRGs, the SDSS QSOs, and the SDSS type 2 quasars are statistically 
different samples in the color-magnitude two-dimensional space. 

• We study the relations among morphology, optical color, and far- infrared luminosity of 
the ULIRGs. These results are broadly consistent with the scenario that the blue optical colors 
originate from more disturbed systems, and these systems are not during the phase when the far- 
infrared emission peaks. However, the morphology parameters are affected by the measurement 
uncertainties, as well as physical properties of the merger. We thus do not conclude any strong 
relations among morphology, far-infrared luminosity, and optical color of these system. 

Selected by their true infrared luminosity with no SED extrapolations, and with the uniformity 
of the SDSS optical photometry, our ULIRG sample is a best low-redshift comparison sample to 
study the color-magnitude relation of high-redshift ULIRGs and sub- millimeter galaxies. Although 
the morphological study of our sample is limited by the physical resolution of SDSS, the size and 
structures of a galaxy might also be smaller at high redshift, and therefore our sample also provides 
a comparison to study the morphology of high-redshift ULIRGs and sub-millimeter galaxies. We 
will discuss the color-magnitude relation and morphology for a sample of the most luminous infrared 
galaxies at z ~ 1 in a succeeding paper. 
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Fig. 1. — Mosaic of RGB color images of all 54 sources. North is up and east is left for each 
image. Each image is 100 pixels across, or ~ 40" at SDSS's pixel scale, and centered on the source. 
The images are arranged according to their g— band G and M20 (see $5]), so that for each column 
G increases towards the top, and for each row M20 increases towards the left. The G and M20 
are printed on each image. A clear trend is seen: more disturbed and multiple nuclei sources are 
located towards the upper left corner of the mosaic. 
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Fig. 2. — Color-magnitude diagram derived from SDSS. Scattered gray dots represent NYU VAGC 
DR4 galaxies with known spectroscopic redshifts within 0.018 < z < 0.265, the same redshift range 
of the ULIRGs. Contours represent 10 equally spaced levels between the minimum and maximum 
number density. Filled circles and stars represent the ULIRG sample, with circles representing 
AGNs and stars representing Hn-like and LINERs, respectively. The solid line is taken from 
Weinmann et al. (2006) as an empirical separation between the red sequence and the blue clouds. 
Source FSC12265 + 0219 and FSC125A0 + 5708 are saturated in the g and r band and open 
circles and arrows are plotted to denote the lower limits. The error bar shows typical photometric 
uncertainties for ULIRGs only, and is dominated by aperture differences among different bands. 
See text for more detailed descriptions on the uncertainties. Typical photometric uncertainties for 
the SDSS comparison sample are less than 0.03 mag in g and r— band (York et al. 2000). The plot 
shows that 24, or ~43% of the ULIRGs lie outside the 90% normal galaxy contour, ~ 46% have the 
typical blue-cloud color, and only ~7% fall in the red sequence. Strikingly, only 3 (~ 6%) ULIRGs 
are located in the green valley, and none of which is an AGN ULIRGs. See text for details. 
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Fig. 3. — K-corrected absolute magnitude (Mo.i r — 51og/i) distributions for (a) the ULIRG sample 
(beige solid line) and the SDSS comparison sample galaxies (black solid line); (b) the AGN ULIRGs 
(solid line) and non-AGN ULIRGs (dashed line); (c) the ULIRG sample (beige dashed line) and 
the three sub-samples of the SDSS comparison sample (solid lines: the red sequence, red; the 
green valley, green; and the blue cloud, blue); (d) the two subsamples of the ULIRG sample (AGN 
ULIRGs: beige solid line; non-AGN ULIRGs: beige dashed line) and the three sub-samples of 
the SDSS comparison sample (the same symbols as in (c)). Note that there is very little overlap 
between the AGN ULIRGs and the green valley in magnitude distribution. 
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Fig. 4. — 01 g — r color distributions for (a) the ULIRG sample and the SDSS comparison sample 
galaxies within the same magnitude range (—19.8 < Mo.i r — 5logh < —23.2); (6) the AGN ULIRGs 
and non-AGN ULIRGs; (c) the ULIRG sample and the three sub-samples of the SDSS comparison 
sample (the red sequence, green valley, and blue cloud) within the same magnitude range (—19.8 < 
Mo.i r — hlogh < —23.2); (d) the two subsamples of the ULIRG sample and the three sub-samples of 
the SDSS comparison sample within the same magnitude range (—19.8 < Mo.i r — 5logh < —23.2). 
The symbols are the same as in Fig. [3j Note that there is little overlap between the AGN ULIRGs 
and the green valley in color distribution. 
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Fig. 5. — Color-magnitude diagram for the 12 unsaturated AGN ULIRGs (filled circles) and the 
residuals after the point source subtraction (empty circles), each pair connecting with a solid line, 
superimposed on the number density contours of the SDSS comparison sample galaxies. The area 
between the two straight lines is the green valley, left: maximum subtraction method; right: smooth 
host subtraction method. See §4.41 for details about the subtraction methods. The subtracted point 
sources contribute on average a small fraction to the total flux, and the residuals are on average 
only 0.45 and 0.47 mag fainter than their un-subtracted counterparts in g and r band, using the 
maximum subtraction method, and 0.27 and 0.31 mag fainter in g and r band, using the smooth 
host method, respectively. Although some sources do appear to be much redder or bluer after the 
point sources being subtracted, the average 01 g — 0,1 r color is almost the same. The residuals are 
barely overlapped with the SDSS comparison galaxies and we do not observe an concentration of 
AGN host galaxies in the green valley. 
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Fig. 6. — Black contours and symbols represent the SDSS field galaxies and the ULIRGs, respec- 
tively. Grey contours and dots represent the SDSS QSOs within the same redshift range of the 
ULIRGs, and blue contours and dots represent the SDSS type 2 quasars within the same redshift 
range. The ULIRGs have comparable optical luminosity to the low-redshift QSOs, but the QSOs 
extend to much bluer regions in the color-magnitude diagram. The ULIRGs are more luminous 
and redder than the SDSS type 2 quasars. Strikingly, the distribution of the SDSS type 2 quasars 
peaks at the green valley. 
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Fig. 7. — Evolutionary tracks are superimposed on the color-magnitude diagram, where the con- 
tours and symbols are the same as in Figj2j Different tracks represent different star formation 
histories: normalized single stellar burst (SSP, solid line; see §4.61 for details), constant star forma- 
tion (dashed line) , star formation histories predicted by numerical simulations of the merger of two 
spiral galaxies (Springel et al. 2005), each galaxy has a dynamical mass of 3.85 x 1O 12 M , with 
(long dashed line) and without (dotted line) AGN activity. The color and magnitude at 1, 2, and 
5 Gyr are marked by stars, triangles, and squares, respectively, for each of the four tracks. The 
evolutionary tracks are calculated using Bruzual & Chariot (2003) with the Salpeter IMF and solar 
metallicity. 
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Fig. 8. — Evolutionary tracks of a combination of a merger-triggered starburst and its passively 
evolved progenitors, with different metallicities and dust extinctions. See £14.61 for more details. 
Three symbols are plotted on each track to mark lGyr, 2 Gyr, and 5Gyr after the starburst. We 
find that dust extinction mainly affects blue optical colors at early ages of the stellar populations, 
and has little effect after 1 Gyr. The metallicity affects both the color range of the stellar populations 
and the evolutionary timescale. The stellar populations with the highest metallicity (bottom panels) 
evolve to the red sequence twice as fast as the lowest metallicity ones (upper panels). 
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Fig. 9.— G - M 20 of the ULIRG sample shown with circles (AGN ULIRGs) and stars (non-AGN 
ULIRGs), superimposed on the comparison sample (dots and contours) in g, r, i, and z bands. 
The comparison sample consists of galaxies on the same images as the ULIRGs, excluding possible 
stellar objects (see text for detailed descriptions). Only sources with signal-to-noise ratio per pixel 
greater than 2.5, and Petrosian radius greater than twice the size of the PSF FWHM (~ 2.6") are 
plotted. The solid line, dashed line and dotted line divide the parameter space into four different 
regions (as marked in the g-band panel) according to LPM04: merger, elliptical, irregular and spiral. 
AGN ULIRGs are not more preferentially concentrated in the merger (non-merger) region than the 
non-AGN ULIRGs. Error bars plotted in the upper left panel show the maximum uncertainties 
in G and M20 estimated from our simulations. The vector in the same panel indicates the upper 
limits of systematic changes in the measurements of G and M20 given increasing background noise 
until the galaxy is no longer detected (either r p < 2 x FWHM, or (S/N) < 2.5). Please refer to 
£|5,3I for a more detailed description. 
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Fig. 10. — Distribution of morphology classifications using G — M20 at each merger stage defined 
in Veilleux et al. (2002) for the same sources. The merging stages are: Tpl - triple mergers, Ilia 
- wide binary pre-mergers, Illb - close binary pre-mergers, IVa - compact mergersm, IVb - diffuse 
mergers, V - old mergers. Values in the horizontal axis stand for the four G — M20 classifications: 
M (merger), E (elliptical), I (irregular), and S (Spiral). 
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Fig. 11. — G-M20 measured in g— band. Different symbol sizes representing different far-infrared 
luminosities, linearly in logarithm scale. Different symbol colors representing different 01 g — 0,1 r 
colors. The small panel shows the distributions of more disturbed systems {G > 0.55) and of less 
disturbed systems (G < 0.55). Sources with redder optical colors are slightly preferentially located 
in less disturbed regions with lower G values. The optically redder sources are slightly less luminous 
in FIR than the bluer ones. All AGN ULIRGs have G > 0.55. However these differences are not 
significant based on the uncertainties of morphological parameter measurements, and the physical 
complexity of the merging systems. 
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Table 1 — Continued 
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Note. — (1): IRAS Faint Source Catalog name (Moshir ct al. 1992). (2): Optical redshifts (Kim & Sanders 1998). (3): FIR luminosity (from 60^m and 100/rni IRAS fluxes; 
Kim & Sanders 1998). (4)-(13): Photometry results, this work (see text). (14): Absolute magnitude in 01 r (see text). (15): 01 g — 0,1 r color (see text). (16): Two Petrosian 
radii in r— band adopted as the photometry aperture (see text), f: For the sources with multiple parts well separated, the apertures for all parts are shown. (17): AGN types by 
previous optical and near-infrared spectroscopic studies. Syl and Sy2 stand for Seyfert 1 and Seyfert 2 types, respectively. References arc: [a] VKS99; [b] Vcron-Cctty & Vcron 
1989; [c] Veilleux et al. 1995; [d] Kim, Veilleux & Sanders 1998; [e] Low et al. 1998. 
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Table 2. Photometry Statistics of the ULIRGs and SDSS Comparison Sample 



ULIRGs 



SDSS Comparison Sample 



Number of Sources 
(Mo.i r ) ± 0M o . lr 
1 9-°- 1 r)±a .i g r _ a . 



Number of Sources 
(M .i r ) ± <t Mo i 
( - 1 g- 01 r)±ao. lg _ . lr 



52™ 
-21.4 ±0.72 
0.57 ±0.22 



436762 
-20.2 ± 1.17 
0.77 ±0.25 



AGN ULIRGs non-AGN ULIRGs Red Sequence Green Valley Blue Cloud 



-22.1 ± 0.60 
0.61 ±0.30 



40 

-21.1 ±0.62 
0.58 ±0.17 



216648 
-20.6 ± 1.01 
0.95 ±0.15 



51904 
-20.4 ± 1.01 
0.76 ± 0.04 



168119 
-19.8 ± 1.26 
0.52 ± 0.18 



I 1 ' The ULIRGs sample excludes the two saturated sources FSC12540±5708 and FSC12265±0219, and these two 
sources arc also AGN ULIRGs. 



Table 3. Two Dimensional Kolmogorov-Smirnov Test for g — r and Mo.i r 





Number 


SDSS QSOs 


Type 2 quasars 




ULIRGs 




AGN ULIRGs 




non- 


AGN ULIRGs 




P(> Z)M 


D 


P(> Z) 


D 


P(> Z) 




D 


P(> Z) 




D 


P(> Z) 


SDSS QSOs 


1170 






0.86 


0.00 


0.59 


1.52 x 10" 


12 


0.79 


2.70 x 10" 


■7 


0.62 


8.31 x 10" 11 


Type 2 quasars 


402 


0.86 


0.00 






0.69 


1.26 x 10" 


IS 


0.87 


1.87 x 10" 


■8 


0.66 


2.11 x 10" 11 


ULIRGs 


54 


0.59 


1.52 x 10~ 12 


0.69 


1.26 x lO" 15 








0.49 


1.17 x 10" 


-2 


0.17 


6.13X10- 1 


AGN ULIRGs 


14 


0.79 


2.70 x 10~ 7 


0.87 


1.87 x 10~ 8 


0.49 


1.17 x 10" 


-2 








0.66 


2.82 x 10~ 4 


non-AGN ULIRGs 


40 


0.62 


8.31 x 10- 11 


0.66 


2.11 x lO" 11 


0.17 


6.13x10" 


1 


0.66 


2.82 x 10" 


-4 






Red Sequence 


216648 


0.98 


0.00 


0.60 


0.0 


0.94 


6.05 x 10- 


32 


0.91 


1.18 x 10" 


-9 


0.94 


2.24 x 10~ 24 


Blue Cloud 


220023 


0.87 


0.00 


0.43 


9.81 x 10~ 45 


0.64 


1.47 x 10- 


IS 


0.80 


1.57 x 10" 


•7 


0.62 


5.62 x 10" 11 



W D is the maximum absolute difference between the two samples' cumulative distribution functions. In two dimensional KS test we considered all 
four possible ranking combinations. For details see Peacock (1983). 

l 2 'Z = J n^~+^ D is the test statistics, where N\ and N2 arc the sample sizes. P(> Z) gives the probability of not rejecting the null hypothesis 
Ho : the two samples are drawn from the same population. The calculation of P(> Z) utilize the function probks in the Numerical Recipes (Press et 
al. 1992). 



